Quantum Melting and Absence of Bose-Einstein Condensation in Two-Dimensional 

Vortex Matter 
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We demonstrate that quantum fluctuations suppress Bose-Einstein condensation of quasi-two- 
dimensional bosons in a rapidly-rotating trap. Our conclusions rest in part on the derivation of 
an exact expression for the boson action in terms of vortex position coordinates, and in part on a 
solution of the weakly-interacting-boson Bogoliubov equations, which simplify in the rapid-rotation 
limit. We obtain analytic expressions for the collective-excitation dispersion, which is quadratic 
rather than linear. Our estimates for the boson filling factor at which the vortex lattice melts are 
consistent with recent exact-diagonalization calculations. 
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vector (g) dependence that is quadratic at small q, in 
contrast to the linear behavior ordinarily characteristic 
of interacting boson systems. The quadratic energy dis- 
persion has fundamental consequences for the nature of 
the interacting-boson ground state: for any v, the inte- 
gral for the fraction of particles outside the condensate 
diverges logarithmically at small q. Thus, Bose-Einstein 
condensation does not occur even at zero temperature in 
the thermodynamic limit. Fluctuations need not cause 
the vortex lattice to melt, however, and our results for 
the melting filling factor are consistent with recent exact- 
diagonalization estimates. 

A system of interacting bosons of mass M in a cylin- 
drical trap (with radial and axial trap frequencies VLr 
and r^z) that is rotating with angular velocity fiz is well- 
described by the rotating-frame Hamiltonian |^,^ 



Recent experiments [Ig] in rapidly-rotating Bose- 
condensed alkali gases [3| have established the occur- 
rence of large vortex arrays. These observations have 
raised a number of fundamental questions about boson 
vortex matter that are suggested by loose analogies with 
the properties of type-II superconductors in a magnetic 
field, and by the lore of the boson quantum Hall effect 
Q. The simplicity of dilute alkali gases, and the arse- 
nal of techniques that have been developed to manipu- 
late them experimentally, make these systems ideal for 
the study of condensed-matter systems containing vor- 
tices. The proportionality between rotation frequency 
and rotating-frame effective magnetic field suggests ||] 
that there will be a maximum rotation frequency, f2c2j 
beyond which vortex-lattice states cannot occur, even 
at zero temperature. This expectation is consistent 
1^ with our understanding from quantum-Hall studies 
that quasi-two-dimensional (2D) charged-boson systems 
in strong magnetic fields form strongly-correlated fluid 
ground states that do not break translational symme- 
try. In theoretical studies of rapidly-rotating boson sys- 
tems, exact-diagonalization studies tend to suggest that 
the bosons form fluid states, while Thomas-Fermi and 
mean-field studies predict vortex-lattice states. Very re- 
cently, it has been suggested 10,^ , on the basis of exact- 
diagonalization studies using periodic boundary condi- 
tions, that the vortex-lattice state of quasi-2D bosons 
melts due to quantum fluctuations when the boson fill- 
ing factor v, the ratio of boson density to vortex density, 
is smaller than ~ 6. 

In this Letter, we present a theoretical study of quan- 
tum fluctuations in rotating Bose-Einstein condensates 
that is based partly on a long-wavelength effective La- 
grangian derived from the microscopic action, and partly 
on microscopic Bogoliubov-approximation calculations in 
the rapid-rotation limit. We find that the collective ex- 
citation energy of quasi-2D rotating bosons has a wave- 
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where V{Y)—g5{v) is a hard-core interaction potential 
of strength g. This Hamiltonian is equivalent to that 
of a system of charge-Q bosons with weakened radial 
confinement, under the influence of a magnetic field 
B = {2MVI/ Q)z. The theory simplifies in the quasi-2D, 
rapid-rotation, weak-interaction limit, in which the in- 
equalities SI ^ -^/Sl^ — > 0, Tig < Q.Z and ng < 2S1 are 
satisfied. (Here n is the boson number density.) All the 
calculations discussed here use the rotating-frame Hamil- 
tonian and assume the rapid-rotation limit in which only 
single-particle states in the lowest Landau level (LLL, 
those states that have the minimum quantized kinetic 
energy) are relevant. (For weakly-interacting bosons, fill- 
ing factors v larger than one are compatible with the LLL 
constraint.) However, as we explain later, we believe that 
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our most important conclusions are more general. In the 
quasi-2D weak- interaction limit, properties of the system 
depend non-trivially only on the boson filling factor 



V = N/Nv = 2Tii^N/A = {N/A){2Trh/2Mn). 
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Here A is the cross-sectional area of the rotating system 
of N bosons, A'^v' is the number of vortices in the system, 
and l—^fL/2MVl is the effective magnetic length. In the 
following paragraphs, we consider unconfined (fi^Slr) 
but interacting 2D bosons in the LLL We discuss 
quantum fluctuations in the vortex-lattice ground state 
by first using a coherent-state path-integral description 
in terms of vortex-position coordinates, which makes the 
physical origin of our surprising findings apparent. We 
then work directly with the boson coordinates by using a 
Bogoliubov approximation, which is a simpler approach 
for quantative analysis of the weak-fiuctuation limit. 

The interacting-boson partition function is Z = 
/ I?[(/)(r)(/)(r)] exp(— where the action is 
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\i5ndr5e + ^{5nfY (3) 



dr i dr 



In the approximate Gaussian-fluctuation form of this ac- 
tion that we employ, the amplitude and phase contri- 
butions are separated: (}> = ^Jne^^ . The usual gradient 
terms are absent in Eq. (||) because of the restriction to 
the LLL, and we have dropped the constant boson kinetic 
energy, Nhfl. The imaginary-time-dependent coherent- 
state label 0(r, r) satisfies the LLL constraint, which 
guarantees that it is completely specified by the locations 
of its zeroes, up to an overall constant. In the mean-field 
ground state, the zeroes of (j> are located at the sites of 
an Abrikosov triangular vortex lattice. We express the 
action in terms of fluctuating vortex positions, u^, by 
writing as a symmetric-gauge LLL wavefunction [Q: 

Nv 

<I>M = *u - ^» - exp(-zz/4£2) = V^e*«, (4) 

where z = x + iy. Here the Zi are the triangular-lattice 
vortex sites, and the vortices of 4>[u] are located at z = 
Zi + Ui'. i.e., Ui is the displacement of the z-th vortex. It 
follows from the LLL constraint that the order-parameter 
phase and amplitude fluctuations are given to leading 
order in u by 
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Fourier-expanding the fluctuation action, we find that at 
long wavelengths it assumes the form 
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where UL = [qxUx{(l, t) -f qyUy{q, T)]/q is the longitudinal 
vortex displacement, UT={qxUy — qyUx)/q is the trans- 
verse vortex displacement, cge is a shear modulus that 
appears here as a phenomenological coefficient (but see 
below), al=2r^(G=0)/^^ a2^2gn'^ {G=Q) / , \Mz)? = 
Eg "(G) exp(zG • r), \Mz)\^ = Eg "'(G) exp(zG • r), 
and G denotes the reciprocal-lattice vectors of the trian- 
gular lattice [0 . 

A closely-related classical energy functional arises in 
the theory of thermal vortex-line fiuctuations in type- 
II superconductors, and leads in that case to the loss 
of long-range coherence p"^-p^ at finite temperatures. 
Here, it follows from our analysis that true off-diagonal 
long-range order, and hence Bose-Einstein condensa- 
tion, is destroyed at zero temperature by quantum fluc- 
tuations for any value of i^, since (|6'(q)p) cx l/q^, 
characteristic of 2D superfluidity with quasi-long-range 
(algebraically-decaying) off-diagonal order. Quantum 
fluctuations do not necessarily melt the vortex lattice, 
however, since, according to the above action, (|u(q)p) 
goes to a constant as (?— >0, impying that the vortex 
lattice can still retain true positional long-range or- 
der. The diverging phase fluctuations of this action 
are accompanied by a gapless collective mode with 
quadratic dispersion, Ecff{q) ~ \/2/3ab5C66('?^)^, where 
/3^B=n2(G=0)/(n(G=0))2=1.1596. To estimate the fill- 
ing factor at which the vortex lattice melts due to quan- 
tum fiuctuations, we need to evaluate the shear con- 
stant cgg, which we now extract from a fully-microscopic 
Bogoliubov-approximation calculation. 

The rotating-frame Hamiltonian [Eq. (^] in second- 
quantized form is 
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where we have switched to the Landau gauge for a 
more convenient microscopic description of vortex-lattice 
states and creates a Landau-gauge periodic- 

boundary-condition single-particle state 
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where Xg — X + sLx with < X < L^, and XLy/i"^ is an 
integer multiple of 27r. The interacting-boson partition 
function may be expressed as a coherent-state path in- 
tegral in which the LLL constraint is imposed explicitly: 
Z = J exp(— S'[i/)]), where 



dr ^ i^xidr - fJ.)ipx +2 E 
{Xi',X2'\V\X^,X2)^x,4x,-i^x,i^xA 
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is the action. Minimizing this action leads to the strong- 
field limit of the Gross-Pitaevskii (GP) equation, 
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A similar equation arises near the upper critical field 
in the Ginzburg-Landau theory of type-II superconduc- 
tors; the role played here by the chemical potential (/i) 
is assumed in that case by the temperature, measured 
from its mean-field-theory critical value. The global min- 
imum of Eq. (H) is known to occur for the triangular 
vortex-lattice state, which has the following form when 
expressed in terms of Landau-gauge single-particle states: 
-f^r) = ExV-^^xW, where = CA[E„^x,„a^ + 
*'5x,(2n+i)aA/2] ) ^^'^ "^A = ill / ^/SP . The Contribution 
to the filling factor from condensate bosons, 1^0=0 
related to /i as described below. 

Expanding the action to second order around the so- 
lution of the GP equation yields the Bogoliubov approx- 
imation for interacting bosons. In the strong-field limit. 
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For the hard-core model considered here, the direct and 
exchange terms in E{X' , X) are identical. The action can 
be further simplified by using the translational symmetry 
of the vortex lattice to transform to a magnetic-Bloch- 
state representation. Then 
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where 6=(gt'o/47ri'^)y 2/V3. The quadratic action in 
Eq. has the familiar form that is diagonalized 
by a Bogoliubov transformation It follows that 

the elementary collective-excitation energy is given by 
Eq = \JS,\ ~ l^qP: ^-iid that the filling factor for bosons 
outside the condensate is 
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+ g^^q^l'^q + V'-qCqV'-q + AqV'qV'-q + AqV'-qV'q) 



where the wave-vector is summed over the triangular- 
lattice Brillouin zone, the prime indicates that q=0 is 
excluded from the sum, and — 2eA(q) ^ fJ-, with 
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We find, in agreement with our long-wavelength anal- 
ysis, that the collective-excitation energy has a de- 
pendence at small q {Eq^0.0836gi^oq'^ for q£^0). As 
a consequence, the density of particles outside the con- 
densate diverges in the thermodynamic limit, consistent 
with divergent condensate phase fluctuations. As in the 
effective-action theory, the vortex-lattice-state density- 
wave order parameter p(G) lacks divergent fluctuations. 
Comparing expressions for the collective-excitation en- 
ergy as q^O yields cee ~ (0.0836)2^1/2/(2/3^^^*) and, 
upon applying a Debye-like model to extend the effective 
theory to shorter wavelengths, ~ It 

then follows from the Lindemann criterion ||l^ that the 
vortex lattice melts due to quantum fluctuations for i> 
smaller than ~ 8, consistent with exact-diagonalization 
estimates 0. 

The rotation frequency at which Bose-Einstein con- 
densation will no longer be evident in realistic finite-size 
trapped-boson systems can be estimated using Eq. (|^). 
In the rapid-rotation limit, (v-vq) is independent of 
the interaction strength but depends on Ny through 
the system-size infrared cutoff, which we estimate as 
9min=l/V27rlVyl2. xhc resulting plot of (I'-iy'o) versus 
system size (i.e., Ny) is shown in Fig. |l|. Adopting 
(v—vo) < uo/2 as a practical criterion for experimentally- 
evident Bose-Einstein condensation leads to the dot- 
dashed line in Fig. ^. The experimental consequences 
of enhanced quantum phase fluctuations in finite-size ro- 
tating traps will be addressed at greater length in future 
work. 

Although our discussion up to this point has been re- 
(6)stricted to the rapid-rotation limit, we believe that our 
main conclusions are more general, because excitations to 
higher Landau levels are massive and become irrelevant 
|]l2| at low energies. The effect of quantum fiuctuations 
on the condensate fraction at moderate rotation rates can 
be roughly estimated by introducing an ultraviolet cutoff 
(7) in the fiuctuation integral in Eq. (^) at Eq^^^ — 2?if7: 
i.e., keeping only fluctuations in which the boson kinetic 
energy does not increase. Calculations of (j^-z^o) for dif- 
ferent values of Eq^^^ are shown in Fig. |l|. In the vortex 
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systems studied in recent experiments [Q, 2|, A'^y^lOO and 
^ 10''— 10^; for these systems, (i/ — i'o)/i'o ^ 10^"*, by 
this crude estimate that neglects fluctuations of the bo- 
son kinetic energy. We conclude that the strengthening 
of quantum fluctuations that we predict for the rapid- 
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rotation limit does not play an important role in these 
experiments. Vortex-fluid states will tend to occur as 
U^^r- i-e., when the rotating- frame trap frequency van- 
ishes. At sufficiently small values of v, we expect these 
vortex-fluid states to Wigner-crystallize into lattices of 
bosons rather than lattices of vortices. Finally, we re- 
mark that confinement in trapped-boson systems will in- 
duce a gap in the collective-excitation spectrum, which 
will reduce the importance of quantum fluctuations. 

In conclusion, we have shown that the zero-point 
motion of the boson vortex lattice in a large rapidly- 
rotating quasi-2D trap eliminates Bose-Einstein conden- 
sation even at zero temperature, in the thermodynamic 
limit. As a result, the condensate fraction in a finite-size 
trapped-atom system may be relatively small, even when 
the interactions between bosons are weak. We predict 
that the fraction of bosons outside the condensate will 
grow logarithmically with the system size, but that the 
vortex-lattice will not necessarily be melted by such fluc- 
tuations. Our estimate for the filling factor below which 
the lattice melts due to quantum fluctuations, i^^S, is 
consistent with exact-diagonalization estimates . Our 
estimate of the size (Ny) at which Bose-Einstein conden- 
sation is lost, shown in Fig. could be tested in current 
traps if the boson density could be reduced to reach the 
low filling factors at which quantum-fiuctuation effects 
become important. Our predictions add to the funda- 
mental interest ^ in trapped-boson vortex-matter stud- 
ies. 
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FIG. 1. Filling-fraction contribution from bosons outside 
the condensate, I'-vo, versus Nv- The effect of Landau- level 
mixing is estimated by including only collective fluctuations 
with energy smaller than the Landau-level separation 2hn. 
The three curves are for 7 = 2£^hQ/{guo) > 0.1 (solid line), 
7=0.01 (dotted line), and 7=0.001 (dashed line). The mini- 
mum condensate filling fraction i/g, below which the conden- 
sate disappears, is shown as the dash-dotted line. 
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